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ABSTRACT 

We present an analysis of Kitt Peak National Observatory and Lowell Observatory observations of 
comet 103P/Hartley 2 obtained from August through December 2010. The results are then compared 
with contemporaneous observations made by the EPOXI spacecraft. Each ground-based dataset has 
previously been investigated individually; the combined dataset has complementary coverage that 
reduces the time between observing runs and allows us to determine additional apparent periods 
at intermediate times. We compare CN coma morphology between ground-based datasets, making 
nine new measurements of apparent periods. The first five are consistent with the roughly linearly 
increasing apparent period during the apparition found by previous authors. The final four suggest 
that the change in apparent period slowed or stopped by late November. We also measure an inner 
coma lightcurve in both CN and R-band ground-based images, finding a single-peaked lightcurve 
which repeats in phase with the coma morphology. The apparent period from the lightcurve had 
significantly larger uncertainties than from the coma morphology, but varied over the apparition in a 
similar manner. Our ground-based lightcurve aligns with the published EPOXI lightcurve, indicating 
that the lightcurve represents changing activity rather than viewing geometry of structures in the 
coma. The EPOXI lightcurve can best be phased by a triple-peaked period near 54-55 hr that 
increases from October to November. This phasing reveals that the spacing between maxima is not 
constant, and that the overall lightcurve shape evolves from one triple-peaked cycle to the next. These 
behaviors suggest that much of the scatter in apparent periods derived from ground-based datasets 
acquired at similar epochs are likely due to limited sampling of the data. 

Subject headings: comets: general — comets: individual (103P/Hartley 2 ) — methods: data analysis 
— methods: observational 


1. INTRODUCTION 

Comet 103P/Hartley 2 was extensively studied during 
the second half of 2010 due to it being the target of the 
EPOXI flyby on 2010 November 4 as well as having a 
close approach to Earth with favorable viewing geom¬ 
etry. EPOXI spacecraft observations revealed that the 
nucleus of Hartley 2 was much smaller than previously 
thought, with its apparent “hyperactivity” du e to a pop¬ 
ulation of icy grains in the coma (e.g., A’Hearn et al 


20111 Hermal yn et al.|20l3l|Kelley et al.|2013| |Frotopapa 

et al. 2014p . Details of the nucle us’s shape and surface 


properties were determined (e.g., |Groussin et al. 2013 


Li et al. 


2013 


Thomas et al. 


2013) and used to explain 


the origins of its jet activity (e.g., |Belton||2013| [Bruck 
Syal et al. 2013). Earth-based observations over mul 


tiple"orbits revealed a secular decrease in activity since 
1991 and a strong asymmetry in production rates around 
perihelion (Combi et al. 2011; Knig htfc Schleicher] 2 013|). 


Many larger-scale properties at or near the time of the 


ocean-like D/H ratio ( 
els and relative abund 

Hartogh et al.||2011), activity lev- 

ances of various gas species (e.g., 

Dello Russo et al.|2011 

Mumma et al. 20111 Weaver et al. 
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2011 ! ICrovisier et al. 2013; |McKay et al. 20131 Boissier 

etah 2014 

), and coma morpology (e.g., |Knight & Schle- 

icher 2011 

Lara et al. 2011; Samarasinha et al. 20 il; 


measurements of the apparent period were determined 
and are discussed in detail below. 

We obtained independent observations of Hartley 2 
at Kitt Peak National Observatory (KPNO) and Low¬ 
ell Observatory on complementary epochs from August 
2010 through January 2011 which have already yielded 
much novel information about Hartley 2 . Data from each 
observatory were analyzed separately, with repetition of 
CN coma features used to determine apparen t periods 


(Knight et al. 2010 Samarasinha et al. 2010). As the 
apparition progressed, continued monitoring of the CN 
coma morphology revealed that the apparent period was 
increasing and suggested that the nucleus was in non¬ 
principal axis (NPA) rotation (Samarasinha et ah]|2010 


2011 Knight & Schleicher 


with coma morp hol ogy 


best matching multiples of three cycles apart (Knight & 
Schleicher || 20 11 ). We later studied the dust morphology, 
identifying two faint sunward-facing dust jets that be¬ 
haved very differently from the strong CN gas jets and 
suggesting that they originated from diffe rent source re¬ 
gions on the nucleus (M ueller et al.||20 13). We analyzed 
the coma morphology ol other gas species, finding that 
the distributions of C 2 and C 3 were similar to that of 
CN, but the distributions of OH and NH were different 
and were likely derived from small ic y grains th at wer e 
subject to radiation pressure (Knight & Schleiche r|2013 ). 
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Numerous observers (A’Hearn et al.| 

20lT| jehin et al. 

2010; 

Knight & Schleicher 2011||lVleech et al. 201T[ bama- 

rasinha et al. 201 If showed that the apparent period near 


usually large amount in only a few months. This change 
cannot always have been linear since the apparent pe¬ 
riod w as nearly unchange d from mid-2009 (|Meech et al. 


2011) until August 2010 (Knight et al. 


201(5]) ; this is un- 


surprising since period changes would be expected to be 
small or nonexistant due to lack of torques w hen activity 
was w eak at large heliocentric d istances. | A’Hearn et al. 
(2011) and Belton et al.| ([2013) identified several peri¬ 
odicities in the EPOXI data ranging from 17-90 hr. In 
addition to a periodicity corresponding to the lengthen¬ 
ing ~17 hr periodicity just discussed, they found a strong 
periodicity near 55 hr (or possibly ^27 hr). This ^55 hr 
periodicity remained constant or increased slightly from 
mid-October until late November, but its behavior in 
September was inconclusive. No evidence for or against 
a change in the longer apparent period has been detected 
via coma morphology. There is no published evidence for 
the extent of any change in the long apparent periods af¬ 
ter the EPOXI observing interval. Understanding the 
behaviors of the short “single-peaked” and long “triple- 
peaked” apparent periods over as wide a range of dates as 
possible is critical for properly interpreting observations 
of Hartley 2 since it may help to reveal the underlying 
component periods of NPA rotatiorj^] 

Despite the wealth of studies on Hartley 2, many 
uncertainties remain, and continued mining of existing 
datasets is worth pursuing. Hartley 2 remains a prime 
target for future spacecraft missions and many of the 
behaviors it exhibited were also seen in another promi¬ 
nent s pacecraft target, IP/Ha lley (cf. |Millis & Schleicher 
1986 Schleicher et al. 1990). It has been w ell e stab¬ 
lished that Hartley 2 is in JNIPA rotation (e.g., A’Hearn 


et al. 2011), but the mode of rotation, long-axis mode 
(LAM) versus short-axis mode (SA M), is not yet defini¬ 


tively known. Belton et al. 


while Samarasinha et al. "(2d 


J2013| argued for a LAM 
Z| pointed out that a SAM 


cannot be excluded. While determining the specifics of 
the NPA rotation are beyond the scope of this paper, our 
ultimate goal is to derive the specific rotational state and 
component periods of NPA rotation. Therefore, in this 
paper we are investigating apparent periodicities in the 
coma morphology (by observing the motion and varia¬ 
tion of morphological features and looking for repetitions 
of patterns) and in the lightcurve (by measuring bright¬ 
ness variations of the inner coma) and how they evolve 
with time in order to constrain the rotational state. The 
apparent repeatability of morphology and/or lightcurve 
features is not necessarily a direct representation of the 
“rotation” of the nucleus, but may be caused by the ori¬ 
entation of the source region(s) on the nucleus appearing 
nearly the same at two separate intervals of time due to 
a combination of NPA rotation, changes in the compo¬ 
nent periods of NPA rotation, changing viewing geom¬ 
etry, and seasonal changes in activity. Thus, we avoid 
using the term “rotation period” when referring to ap- 


5 Note that there are three component periods of NPA rota¬ 
tion: P. 0 , Pfl, P<£, that P^ and Pq are coupled to each other, and 
that we are using the “L-conve ntion” for defining the Euler angles 
( jSamaras inha & Muelle r|2015| and references therein). 


parent periods. 

In this paper we combine KPNO and Lowell imaging 
data to gain new insight into the rotation of Hartley 2. 
The datasets were temporally offset from one another 
except at the time of the EPOXI encounter, and the 
combined dataset allows us to determine apparent pe¬ 
riods from morphology at intermediate times that are 
impossible to derive for either individual dataset. The 
combined dataset provides additional lightcurve data at 
times beyond the EPOXI lightcurve coverage of Hart¬ 
ley 2 and which have not been reported in the liter¬ 
ature by any other observers. Previous studies have 
yielded systematic offsets between apparent periods de¬ 
termined from lightcurve and coma morphology, e.g., be- 


coma morphology (Jehin et al.|2010 

Knight & Schleicher 

2011 Samarasinha et al. 2011) are higher than those 

from light curves at similar epochs ( 

A’Hearn et al. 2011 

Drahus et al. 2011[ Waniak et al. 

2012; Belton et al. 


of apparent periods for lightcurve and coma morphol¬ 
ogy from the same data to test for differences. We also 
utilize the EPOXI lightcurve data to better understand 
our data and explore possible biases in apparent peri¬ 
ods determined from more sparsely sampled Earth-based 
datasets. 

In Section [2] we briefly summarize the observations 
and reductions of each dataset. In Sections [ 3 ] and 0 ] we 
present studies of the coma morphology and inner coma 
light curves, respectively. Finally, we discuss the implica¬ 
tions of these new results in Section 0 

2. OBSERVATIONS AND REDUCTIONS 


As describ ed previously in Samarasinha et al. (2011) 
and Mueller et al. (2013), observations with the KPJN 0 
2.1-m telescope haa "a field of view of 6.5arcmin x 
10 arcmin and a binned pixel scale of 0.305 arcsec. Lowell 
observations were acquired with the Hall 42-in (1.1 m) 
and 31-in (0.8 m) telescopes, which had square fields of 
view of 25.3 arcmin and 15.7 arcmin, respectively, and 
binned pixel scales of 1.48 arcsec and 0.92 arcsec, respec¬ 
tively (Knight & Schleicher] 2011,2013). All images were 
trailed at the comet's rate of motion. Broadband R and 
HB narrowband comet filters (Farnham et al.|2000) were 
used on all nights, however, we utilize only broadband R 
or narrowband CN images throughout this paper. Im¬ 
ages in additional filters were acquired but generally had 
worse signal to noise than R or CN and were taken with 
insufficient frequency for the studies herein. All 31-in 
data were snapshots and were, therefore, only used for 
coma morphology assessment. Images were processed us¬ 
ing standard bias subtraction and flat fielding techniques. 
The observing circumstances and weather conditions are 
given in Table [l] 

For morphological analysis, we employed enhancement 
techniques to remove the bulk background and make the 
fainter features more obvious. KPNO images were en¬ 
hanced by dividing by each image’s azimuthal average 
profile (with high pixel rejection) and Lowell images were 
enhanced by dividing by each image’s azimuthal median 
profile. These enhancement techniques produce nearly 
identical results, a fact which was confirmed by compar¬ 
ison of enhanced images for all nights on which simulta¬ 
neous data were obtained. We created subimages of the 
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Table 1 

Summary of Hartley 2 imaging observations and geometric parameters in 2010. a 


UT Date 

UT Range 

Tel b 

AT C 

(days) 

r H d 

(AU) 

A e 

(AU) 

e ( 

(°) 

PA Suns 
(°) 

Conditions 

Used in 
Lightcurves? 

Aug 13 

03:28-11:48 

42in 

-75.9 

1.463 

0.567 

30 

22 

Photometric 

Yes 

Aug 14 

03:03-11:48 

42in 

-74.9 

1.454 

0.557 

30 

21 

Photometric 

Yes 

Aug 15 

03:15-12:10 

42in 

-73.9 

1.446 

0.547 

30 

20 

Clouds 

Yes 

Aug 16 

07:40-11:54 

42in 

-72.8 

1.437 

0.536 

31 

19 

Clouds 

Yes 

Aug 17 

07:48-11:01 

42in 

-71.9 

1.421 

0.527 

31 

18 

Clouds 

Yes 

Sep 1 

03:24-11:47 

KP2.1 

-56.9 

1.310 

0.391 

34 

3 

Photometric 

Yes 

Sep 2 

03:22-12:01 

KP2.1 

-55.9 

1.302 

0.383 

34 

3 

Cirrus 

Yes 

Sep 3 

03:08-11:04 

KP2.1 

-55.0 

1.295 

0.375 

35 

2 

Cirrus 

Yes 

Sep 9 

02:53-12:16 

42in 

-48.9 

1.252 

0.328 

37 

357 

Photometric 

Yes 

Sep 10 

02:36-12:09 

42in 

-48.0 

1.245 

0.320 

37 

357 

Photometric 

Yes 

Sep 11 

02:33-12:08 

42in 

-47.0 

1.238 

0.313 

37 

356 

Photometric 

Yes 

Sep 12 

02:33-12:11 

42in 

-46.0 

1.231 

0.305 

38 

356 

Clouds 

Yes 

Sep 13 

02:35-12:11 

42in 

-44.9 

1.224 

0.298 

38 

355 

Clouds 

Yes 

Sep 30 

03:19-12:12 

KP2.1 

-27.9 

1.127 

0.188 

44 

3 

Patchy clouds 

Yes 

Oct 1 

02:25-11:30 

KP2.1 

-27.0 

1.123 

0.183 

45 

5 

Patchy clouds 

Yes 

Oct 2 

03:53-12:24 

KP2.1 

-25.9 

1.118 

0.177 

45 

7 

Patchy clouds 

Yes 

Oct 3 

05:43-12:36 

KP2.1 

-24.9 

1.114 

0.172 

46 

10 

Cloudy 

Yes 

Oct 4 

03:06-09:02 

KP2.1 

-24.0 

1.110 

0.167 

46 

12 

Cirrus 

Yes 

Oct 12 

03:12-12:34 

31in 

-15.9 

1.082 

0.134 

49 

43 

Photometric 

No 

Oct 13 

03:13-12:43 

31in 

-14.9 

1.079 

0.131 

49 

48 

Photometric 

No 

Oct 14 

03:13-12:40 

31in 

-13.9 

1.076 

0.129 

50 

52 

Photometric 

No 

Oct 15 

03:06-05:33 

31in 

-13.1 

1.074 

0.127 

50 

56 

Clouds 

No 

Oct 16 

05:01-12:21 

42in 

-11.9 

1.072 

0.125 

51 

61 

Thin cirrus 

Yes 

Oct 17 

05:00-12:38 

42in 

-10.9 

1.070 

0.123 

51 

65 

Clouds 

Yes 

Oct 19 

10:56-12:24 

42in 

-8.8 

1.066 

0.121 

52 

73 

Clouds 

No 

Oct 31 

07:10-12:36 

31in 

+3.2 

1.060 

0.140 

58 

99 

Thin cirrus 

No 

Nov 1 

07:15-12:45 

31in 

+4.2 

1.060 

0.144 

59 

100 

Thin cirrus 

No 

Nov 2 

07:37-12:35 

KP2.1 

+5.2 

1.061 

0.148 

59 

102 

Photometric 

Yes 

Nov 2 

06:45-12:54 

42in 

+5.2 

1.061 

0.147 

59 

102 

Photometric 

Yes 

Nov 2 

07:27-10:32 

31in 

+5.1 

1.061 

0.147 

59 

102 

Photometric 

No 

Nov 3 

07:47-12:53 

KP2.1 

+6.2 

1.062 

0.152 

59 

103 

Cirrus 

Yes 

Nov 3 

06:41-13:01 

42in 

+6.2 

1.062 

0.151 

59 

103 

Photometric 

Yes 

Nov 4 

07:48-12:49 

KP2.1 

+7.2 

1.063 

0.156 

59 

104 

Cirrus 

Yes 

Nov 4 

06:39-13:07 

42in 

+7.2 

1.063 

0.155 

59 

104 

Thin cirrus 

Yes 

Nov 5 

07:50-12:54 

KP2.1 

+8.2 

1.065 

0.160 

59 

105 

Cirrus 

Yes 

Nov 5 

07:44-10:45 

31in 

+8.1 

1.065 

0.159 

59 

105 

Photometric 

No 

Nov 6 

07:50-12:54 

KP2.1 

+9.2 

1.066 

0.164 

59 

106 

Cirrus 

Yes 

Nov 6 

07:41-10:58 

31in 

+9.1 

1.066 

0.163 

59 

106 

Thin cirrus 

No 

Nov 7 

07:45-12:53 

KP2.1 

+10.2 

1.068 

0.168 

59 

107 

Patchy clouds 

Yes 

Nov 7 

06:48-13:09 

42in 

+10.2 

1.068 

0.168 

59 

107 

Photometric 

Yes 

Nov 8 

08:04-12:46 

KP2.1 

+11.2 

1.070 

0.173 

59 

108 

Photometric 

Yes 

Nov 10 

07:59-13:11 

31in 

+13.2 

1.074 

0.181 

58 

109 

Clouds 

No 

Nov 12 

07:56-12:30 

31in 

+15.2 

1.079 

0.191 

58 

111 

Clouds 

No 

Nov 13 

08:09-13:10 

31in 

+16.2 

1.082 

0.196 

57 

112 

Clouds 

No 

Nov 16 

08:08-13:08 

31in 

+19.2 

1.091 

0.211 

56 

114 

Photometric 

No 

Nov 26 

07:45-12:43 

31in 

+29.2 

1.132 

0.262 

51 

123 

Photometric 

No 

Nov 27 

07:48-12:46 

31in 

+30.2 

1.137 

0.267 

50 

124 

Photometric 

No 

Dec 9 

07:02-13:19 

42in 

+42.2 

1.205 

0.330 

42 

136 

Thin cirrus 

Yes 

Dec 10 

06:48-08:50 

42in 

+43.1 

1.212 

0.336 

41 

137 

Thin cirrus 

Yes 

Dec 11 

07:23-12:11 

KP2.1 

+44.1 

1.218 

0.342 

41 

138 

Cloudy 

Yes 

Dec 12 

07:11-11:59 

KP2.1 

+45.1 

1.225 

0.347 

40 

139 

Cirrus 

Yes 

Dec 13 

07:17-12:02 

KP2.1 

+46.1 

1.231 

0.353 

39 

141 

Cirrus 

Yes 

Dec 14 

06:54-11:58 

KP2.1 

+46.9 

1.238 

0.358 

39 

142 

Photometric 

Yes 

Dec 15 

07:00-11:58 

KP2.1 

+48.1 

1.245 

0.364 

38 

143 

Photometric 

Yes 

Dec 15 

08:57-09:13 

42in 

+48.1 

1.245 

0.363 

38 

143 

Clouds 

No 


a All parameters are given for the midpoint of each night’s observations. 

k KP2.1 = KPNO 2.1-m Telescope, 42in = Lowell Observatory Hall 42-in Telescope (1.1 m), 31in = Lowell Observatory 31-in Telescope (0.8 m). 


c Time relative to perihelion. 

^ Heliocentric distance. 
e Geocentric distance. 
f Solar phase angle. 

g Position angle of the Sun, measured from north through east (counterclockwise when displayed with north up and east to the left). 


enhanced images that were centered on the nucleus and 
had the same physical scale throughout the apparition, 
^50,000 km at the comet. Because of the smaller field of 
view of the KPNO data, the outer portion of this region 
was not always observed and these regions were set to 0. 

Aperture photometry was conducted using two aper¬ 
tures that were fixed in physical size at the comet (radii 
of 1000 and 2000 km) and five apertures that were fixed 
in angular scale (radii of 3, 6, 9, 12, and 15 arcsec). 
Larger apertures had more signal but also diluted the in¬ 
trinsic variability and suffered more contamination from 


background stars. Conversely, smaller apertures reduced 
stellar contamination and showed larger amplitude vari¬ 
ations, but had less signal and, if too small, did not 
necessarily encompass all of the feature(s) causing the 
brightness variation. Ultimately the 9 arcsec radius aper¬ 
ture proved to be the optimum size and was used for the 
light curve analysis. Appropriately large background an¬ 
nuli were utilized for sky subtraction for each run, but 
there was likely always some coma present in the back¬ 
ground annulus. Coma contamination diminishes the 
lightcurve amplitude and affects night to night normal- 
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ization, but does not affect the timing of extrema and 
therefore was of minimal concern. 

We did relative p hotometr y with on-chip stars as out¬ 
lined in Knight et al. (2011) for all Lowell images and 
for KPNO images from December, allowing us to recover 
a useful lightcurve on non-photometric nights for extinc¬ 
tion <0.5 mag. Due to the combination of a small field 
of view and the comet’s high rate of motion, this was 
not possible for KPNO images from September, October, 
and November. Therefore non-photometric nights dur¬ 
ing these months could not be used for R-band lightcurve 
analysis since the uncertainties due to clouds were of 
order the R-band amplitude. We have included CN 
lightcurves from these nights because they have signif¬ 
icantly larger amplitudes, but we were cautious of over¬ 
interpreting unusual lightcurve behavior on these nights. 

3. CN COMA MORPHOLOGY 


As discussed in |S amarasinha et al. (2011) and Knight 
& Schleicher| (|201l|), large scale features were evident m 
the CN coma that changed smoothly during a night and 
varied from night to night. The repetition of these fea¬ 
tures on different nights within a given run previously 
allowed us to determine the apparent period. The CN 
images were not significantly contaminated by either con¬ 
tinuum or other gas species, allowing us to use data from 
non-photometric nights and greatly increasing the num¬ 
ber of nights available for morphological comparison than 
if we only used photometric nights. The only large scale 
feature seen in the R images was the dust tail which var¬ 
ied little from night to night. Two relatively short, sun¬ 
ward facing features were also seen in the highest signal 
to noise images (e.g., Mueller et al. 2013). While their 
orientations changed over diurnal time scales, these vari¬ 
ations were insufficient to constrain the apparent period 
from our two datasets (Samarasinha et al., in prep., are 
attempting to do this using o bservations fro m loca tions 
around the world; see also Samarasinha et al.| 2012) and 
so R-band morphology is not discussed further. 

As with our previous investigations of the apparent pe¬ 
riod within an individual run, we compared the shape, 
size, radial distance, orientation, and relative intensity 
of CN coma features in order to find pairs of images 
with similar morphology. Here we concerned ourselves 
only with images in previously uncompared datasets, e.g. 
Lowell data from August 13-17 with KPNO data from 
September 1-3, KPNO data from September 1-3 with 
Lowell data from September 10-13, etc. We compared 
all images in consecutive datasets and noted all pairs 
of images with similar morphology. Selected pairs are 
shown in Figure^ Given the large baseline between ob¬ 
serving runs and the reasonably high cadence, typically 
^20 min for KPNO and ~60 min for Lowell, interpola¬ 
tion between images was not necessary. 

We determined the apparent period between pairs of 
images by dividing the time between images by the inte¬ 
ger number of cycles that yielded a period closest to the 
expected average period based on extrapolation from pre¬ 
viously published results. Solutions of N+l cycles and 
N—l cycles were also considered for an initial solution 
of N cycles. However, in most cases the 7V±1 cycle solu¬ 
tions had larger standard deviations, yielded periods that 
were beyond the error bars of previously published solu¬ 
tions, and were inconsistent with a smoothly changing 




100813 6:21 42in 100901 11:07 2.1m 



100903 5:36 2.1m 100910 9:19 



Figure 1. Pairs of matching CN images from different observing 
runs. The date (YYMMDD), UT time (HH:MM), and telescope 
are given below each image. The time sequence moves from the left 
to right pair and then from top to bottom. Each image is centered 
on the nucleus, has been enhanced as described in the text, has 
north up and east to the left, and is ~50,000 km across at the 
comet. The grayscale stretch is the same for both images in a pair, 
but is different from pair to pair. White is bright and black is 
dark, and trailed stars appear as diagonal streaks in some images. 
Some 2.1-m images have a circular black border; these are regions 
where no data were obtained due to the smaller field of view. Slight 
mismatches in orientation are primarily due to geometry changes 
between runs, e.g., between October 16 and November 6. 


period. These were therefore considered unlikely. Fur¬ 
thermore, sometimes the N± 1 solutions could be ruled 
out within a single run, e.g., when data were obtained 2- 
3 days apart. The final solution for the period between 
two observing runs was the mean period of all measured 
pairs between the two runs and the 1-cr uncertainty was 
the standard deviation of the measured periods. The 
time at which this mean period applied was the mean 
of the midpoints of all pairs. We present a summary of 
our period solutions in Table [2j including the AT+1 and 
N—l solutions, and plot our preferred solutions along 
with other published periods in Figure [2] The periods 
centered on November 9.40, November 17.45, and De¬ 
cember 4.93 were calculated from three or fewer pairs of 
images. While the formal uncertainty in these periods 
is less than 0.1 hr, we consider the real error to be no 
better than 0.5 hr, and round to the neareast 0.1 hr. 

It was more difficult to match features in the coma 
between observing runs than within individual observ¬ 
ing runs because the longer baseline resulted in changes 
in the ratios of component periods, larger differences in 
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Table 2 

Apparent periods measured from coma morphology 


Date 

Range #1 

Date 

Range #2 

Average 

Date 

Intervening 

Cycles 

Num. 

Pairs 

Period (hr) 

N cycles 

Period (hr) 
N+l cycles 

Period (hr) 
N— 1 cycles 

Aug 13-17 

Sep 

1-3 

Aug 23.87 

22-27 

20 

17.03 

± 

0.03 

16.41 

± 

0.03 

17.70 

± 

0.08 

Sep 1-3 

Sep 

10-13 

Sep 7.04 

10-17 

70 

17.23 

± 

0.08 

15.97 

± 

0.19 

18.71 

± 

0.26 

Sep 10-13 

Sep 

30-Oct 4 

Sep 22.04 

23-32 

36 

17.41 

± 

0.06 

16.82 

± 

0.06 

18.05 

± 

0.10 

Sep 30-Oct 4 

Oct 

12-16 

Oct 8.13 

12-20 

16 

17.97 

± 

0.21 

16.90 

± 

0.35 

19.20 

± 

0.16 

Oct 16-19 

Nov 

6-8 

Oct 28.17 

26-29 

14 

18.38 

± 

0.22 

17.74 

± 

0.19 

19.08 

± 

0.25 

Nov 3-6 

Nov 

10-16 

Nov 9.40 

7-14 

3 a 

18.5 

± 

0.5 

16.6 

± 

0.5 

20.8 

± 

0.5 

Nov 7 

Nov 

27 

Nov 17.45 

26 

2 a 

18.4 

± 

0.5 

17.8 

± 

0.5 

19.2 

± 

0.5 

Nov 27 

Dec 

11-13 

Dec 4.93 

18-21 

2 a 

18.5 

± 

0.5 

17.6 

± 

0.5 

19.5 

± 

0.5 

Dec 9-10 

Dec 

10-13 

Dec 11.24 

1-5 

6 

18.45 

± 

0.23 

14.06 

± 

2.45 

24.05 

± 

0.90 


a Periods given with 0.1 hr accuracy and uncertainty set to 0.5 hr rather than the formal uncertainty (<0.1) to avoid implying misleading precision 
due to small numbers of matched pairs of images. 
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Figure 2. Apparent periods of Hartley 2 as a function of time 
relative to perihelion (AT). Symbols are defined in the legend. 
Coma morpholo gy was used to dete r mine apparent periods in thi s 
wor k (Table [2|>, [J ehin et al.| ( |2010| > , |Knight fc Schleicher| ( |2011| > , 
and |Samarasmha et al.|(|2011|) . App arent periods w e re det e rmined 
from p hotometric measuremen ts by |Drahus et al.| (|2011|), [M eech 
e^aljj^Oll^j A’ Hearn et al.| |2011|> 7 TWaniair et al.| ( |2012| >, | hTcT^| 
ton et al.|(|2013| , and this work (Tabl ej3[ both CN a nd R results 
are plotted with the same symbols). Harmon et al. I poll) used 
radar observations of the nucleus. Error bars are shown but do not 
necessarily have the same significance between different data sets 
due to differing methodologies. Note that the upper error bar on 
the lightcurve point at AT = +44 days is truncated and actually 
extends to 20.5 hr. 


viewing geometry, etc. The coma morphology from mid- 
August through mid-September was reasonably similar. 
This allowed a large number of matching pairs of images 
with a high degree of confidence in the matches, result¬ 
ing in a small standard deviation. However, the coma 
morphology from late September onward (perigee was 
on October 20) varied substantially and resulted in far 
fewer matching features from run to run and significantly 
less certainty about the mean period. Furthermore, as 
the coma features changed from a nearly face-on spiral in 
August to a side-on hour-glass shape by November, in¬ 
creasingly more subsets of images had little to no easily 
distinguishable coma features. The absence of distinct 
coma features meant that such images could not be used 
to determine matching pairs, although they were useful 
for eliminating alternate apparent periods. 

As seen in Figure [2j our newly measured apparent peri¬ 
ods agree well with those previously published. The scat¬ 
ter in our periods compared to the ensemble is smaller 
than it was for either of our individual datasets, presum¬ 
ably due to the longer baseline between image pairs. The 
one exception is that the period determined for late Au¬ 
gust appears somewhat higher than expected relative to 
the individual periods in mid-August and early Septem¬ 


ber, but is likely more robust because of the longer base¬ 
line. These new measurements are also more consistent 
with previous photometry-based periods, suggesting that 
there is not a systematic difference between techniques. 
The post -EPOXI flyby measurements have large uncer¬ 
tainties but, due to the number of epochs, give a stronger 
indication of a likely flattening of the rotation period 
than any previous results. We will return to these topics 
in more detail in Section [U 


4. LIGHTCURVES 

Since Hartley 2 has a small nucleus and was highly 
active throughout our observations, its nucleus signal was 
swamped by the coma signal. The strong CN jets made 
us optimistic that there would be enough variation in 
the CN signal to measure a lightcurve. Furthermore, 
our observing cadences were high since we were looking 
for rotational variation in the CN coma morphology, so 
our CN temporal coverage was good. Since the dust is 
entrained in escaping gas, the potential large variation of 
the CN gas suggested that a corresponding variation in 
the R-band brightness was possible. The R-band images 
also had high signal to noise and good temporal coverage. 

An inner coma lightcurve is typically much more com¬ 
plicated than a nucleus lightcurve, as it may depend on 
factors including the number and relative strength of the 
jets, the location of the jets, the duration of activity 
for each jet due to its local day/night, and pr ojectio n 


effects (e.g. , the lightcurve of lP/Halley by Schleicher 


et al. 1990). Further complicating matters, NPA ro 


tatlon means a jet will not necessarily be in the same 
orientation each “rotational” cycle. EPOXI obtained 
a coma lightcurve in the clear filter from September 5 
through November 26 which exhibited a single-peaked 
sinusoi d that varied in amplitude over three “rotation” 
cycles dA’Hearn et al.|2011||Belton et al.|2013|). Our data 
extend before and after the EPOXI data, but the EPOXI 
findings gave us reason to believe that the lightcurve dur¬ 
ing our observations would exhibit a similar behavior. 

The inner coma photometry is plotted in Figures [ 3 ] and 
[4] Figure [3] shows the lightcurve as a function of UT on 
a night by night basis for each observing run, with R 
and CN plotted on the same panels. Figure [4] shows the 
same data but phased to the best “period” for that run 
or pair of ad jacent runs (the phasing will be discussed in 
Section 4.2). For clarity, the CN and R lightcurves are 
plotted separately in Figure [4j 


4.1. Lightcurve Shape and Amplitude 
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Figure 4. Rotationally phased lightcurves for each observing run or pair of runs. The left column is the CN lightcurve and the right 
column is the R lightcurve. The period to which the data are phased is given in the bottom right of each panel. Zero phase was set to 
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All magnitudes have been normalized to A = 1 AU by an inverse square law. The magnitude range is different from panel to panel, but 
the span is 1.5 dex for all CN panels and 0.6 dex for all R panels. 


As expected, the lightcurve displayed a single-peaked 
shape corresponding to fluctuations in brightness of the 
inner coma. The CN lightcurve displayed a much larger 
amplitude than the R lightcurve. It is difficult to mea¬ 
sure the peak to trough variation precisely because no 
single night included both the lightcurve maximum and 
minimum and because the amplitude appears to vary 
from cycle to cycle. However, we estimate the peak to 
trough variation in a 9 arcsec aperture was at least 0.4- 
0.5 mag throughout the apparition for CN, while for R 
it was less than 0.1 mag in August and September and 
less than 0.2 mag in October, November, and Decem¬ 


ber. These variations are smaller than those reported by 

pH 

mt< 


A’Hearn et al. (2011) for grains, CO 2 , and H 2 O because 
the EPOX1 lightcurves were measured in a much smaller 
aperture (191 km across) that responds to changes over 
very short timescales and results in larger amplitudes. 

The lower amplitude in R relative to CN is not sur¬ 
prising. Gas molecules (e.g., CN) are typically travel¬ 
ing at significantly higher velocities than are dust grains 
(e.g., the majority of the signal in the R images) and 
will therefore leave a given aperture more quickly. Fur¬ 
thermore, dust grains vary in size, and larger grains will 
travel slower than small grains, with the largest grains 
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remaining near the nucleus for a long time. The higher 
velocity of gas and the dispersion of dust velocities com¬ 
bine to reduce the variability of the dust signal in the 
inner coma as compared to the gas. 


was apparently increasing throughout the apparition so 
it was difficult to use a longer baseline to improve the 
apparent period. We will discuss this last point in more 
detail below. 


4.2. Apparent Period 

We determined the apparent period of our lightcurve 
data by plotting all data from a run, or a pair of runs 
if close together, phased to a particular period with zero 
phase set to the time of perihelion (2010 October 28.257), 
then changing the period slightly and redisplaying the 
phased lightcurve in real time. We looked for alignment 
of lightcurve features on different cycles but, owing to the 
varying amplitude and instrinsic brightness from night to 
night even when photometric, extrema needed only agree 
in phase to be considered aligned. We estimated the vi¬ 
able range of apparent periods as the times when the 
extrema first became noticeably out of phase. Visually 
phasing our CN lightcurves yielded periods at five epochs 
between mid-August and mid-December (some subsets 
of data used for coma periods were inconclusive for con¬ 
straining the lightcurve period). Because there was much 
less variation in the R lightcurve, we were only able to 
constrain the period with it at two epochs, in November 
and December. The apparent CN and R lightcurve pe¬ 
riods are given in Table [3] and plotted in Figure [2j the 
data are shown phased to our preferred apparent period 
for each run in Figure [4j We attempted to determine pe¬ 
riods using Fourier analysis, but could not do so reliably 
due to the limitations of the dataset. Our phase cover¬ 
age was insufficient to significantly constrain the longer 
(triple-peaked) period. 

These periods are generally consistent with the peri¬ 
ods we derived from the coma morphology (Section [3]), 
but with larger uncertainties. The larger uncertainty is 
due to a combination of factors. First, we had numer¬ 
ous non-photometric nights which reduced the number 
of potential matches. These nights were generally sal¬ 
vageable by using comparison stars in the field of view. 
However, portions of nights with extinction greater than 
0.5 mag and KPNO nights in October and November 
that had more than minimal cirrus (background stars 
moved quickly across the chip in these months) were 
effectively unusable for constraining the period. Also, 
some increased level of uncertainty was due to inter¬ 
mittently losing individual images because of contami¬ 
nation from background stars in the photometric aper¬ 
ture. More problematic were the quirks of the lightcurve 
itself. The relatively simple single-peaked curve with a 
period near 16-19 hr meant that we could only sample 
a portion of the lightcurve on a given night, and the 
amount sampled decreased during the apparition as the 
nightly observing window shrunk towards the end of the 
year. If that section did not include a clear minimum 
or maximum it was very difficult to use to refine the 
period, although it was useful for excluding some peri¬ 
ods. Furthermore, the lightcurve itself varied in ampli¬ 
tude and intrinsic brightness from one cycle to the next so 
sections without extrema could not be definitively com¬ 
pared. Thus, the lightcurve required many nights of data 
to be sampled sufficiently for phasing. Compounding this 
problem, the shape of lightcurve features was apparently 
not constant fro m on e cycle to the next, as will be dis¬ 
cussed in Section |T4j Finally, the time between extrema 


4.3. Comparison with EPOXI Lightcurve 

Calibrated aperture photometry from EPOXI s 
Medium Resolution Visible CCD (MRI) images is 
available through the Planetary Data System for Oc¬ 
tober 1-November 26, the inte rval when Hartley 2 w as 
observed nearly continuously (Williams et al. 2013 


Thus, we were able to directly compare our inner coma 
lightcurve with that obtained by the spacecraft during 
our October and November runs. EPOXI MRI photom¬ 
etry is available in the OH, ultraviolet continuum, CN, 


C 2 , green continuum, and clear filters (see Hampton 
et al. 2005 for detailed filter information), but only the 
clear filter was sufficiently well sampled for direct com¬ 
parison with our data. Since the EPOXI CN lightcurves 
were undersampled and our R-band lightcurves showed 
much less variation than our CN lightcurves, we chose 
to compare our CN lightcurves with the MRI clear 
filter lightcurve (Figure [ 5 ]). Both ground-based and 
EPOXI data are plottea with a circular aperture of 
radius 6 arcsec. While smaller than the aperture we 
employed for earlier studies, this aperture size was the 
best match to the ensemble of EPOXI data, where 
the spacecraft-centric distance varied from 0.24 AU on 
October 1 (^30% farther than the geocentric distance 
at that time) to <1000 km on November 4. Times 
were corrected for the light travel time and fluxes were 
multiplied by the observer-centric distance (A). The 
flux correction accounts for the changing observer-comet 
distance (a factor of A 2 ) and a conversion from a fixed 
aperture in angular size to a fixed aperture in linear 
size (A -1 , under the assumption that the coma is in 
free expansion so the surface brightness falls off as 
the inverse of the distance to the comet, p _1 , at these 
scales). Due to differences in the bandpasses, it was 
necessary to scale up the EPOXI data by a factor of 6 
to improve their visibility, but this had no effect on the 
timing of features in the lightcurve. 

The ground-based and EPOXI lightcurves behave sim¬ 
ilarly, with segments acquired at the same times exhibit¬ 
ing the same general trends, and gaps in the EPOXI 
dataset being filled well by the ground-based dataset. 
While the ground-based coverage is considerably sparser 
than the spacecraft coverage, the extrema appear to line 
up well in October. There are small systematic shifts be¬ 
tween the datasets in November, with the ground-based 
lightcurve extrema occuring slightly later in November 
due to the significantly larger aperture size at the comet 
(cf. Li et al.||2012|). The November extrema are aligned 
when we plot apertures of the same size at the comet; 
we do not show this because the EPOXI -comet distance 
changed too rapidly during this interval and the plot 
becomes difficult to read. Since the EPOXI spacecraft 
viewed Hartley 2 from a different direction than the 
Earth, the agreement of the shape of the lightcurve and 
the timing of the extrema suggest that the lightcurve 


6 While Belton et al. (2013| analyzed the EPOXI lightcurve 
beginning on September 0 , the earliest data available in the archive 
is from October 1. 
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Table 3 

Apparent periods measured from phased lightcurves 


Date 

Range #1 


Date 

Range #2 


Midpoint 


CN 


R 


Best (hr) 


Range (hr) 


Best (hr) 


Range (hr) 


Aug 13-17 a 
Sep 1-3 
Sep 30-Oct 4 
Nov 2-8 a 
Dec 9-15 a 


Sep 10-13 
Oct 16-17 


Aug 15 
Sep 7 
Oct 8 
Nov 5 
Dec 12 


16.5 
17.0 
18.0 
19.2 

18.5 


16.1- 16.7 
16.9-17.2 
18.0-18.1 
19.0-19.4 

18.1- 18.6 


18.8 

19.0 


18.6-18.9 

17.5-20.5 


Apparent periods determined from a single set of consecutive nights. 



Oct 4 


Oct 6 




Nov 2 


Nov 4 


Nov 6 


Nov 8 


0.0 0.2 


0.4 0.6 

Phase 


0.8 1.0 


Figure 5. KPNO (blue circles) and Lowell (purple triangles) CN 
lightcurves overpl otted wit h the EPOXI (red crosses) clear filter 
light curve from Williams et al. (2013]). EPOXI CN filter photom¬ 
etry is available but was acquired much less frequently than the 
clear filter so we have elected to display the clear filter data. All 
fluxes were extracted in apertures of radius 6 arcsec and were mul¬ 
tiplied by the observer-centric distance (A) to remove the effect of 
the changing geocentric distance. All times were corrected for light 
travel. The EPOXI data were scaled up by 6 X to improve visibility; 
this had no effect on the timing of features in the lightcurve. The 
significantly smaller scale at the comet for EPOXI in November 
results in a much larger amplitude due to the decreased crossing 
time for grains in the aperture and the decrease in smearing out of 
the signal due to velocity disperson. 

shape is due to changes in the production rate. If, in¬ 
stead, the lightcurve shape was due to changing viewing 
geometry of structures in the coma, we would expect an 
offset in the timing of the lightcurve extrema or possibly 
different trends. 

4.4. Phasing of the EPOXI Lightcurve 

Having shown that our ground-based lightcurves 
match the EPOXI lightcurve, we now briefly consider 
the EPOXI lightcurve by itself. We do not intend for 
this to be a thorough analysis of the EPOXI data, but a 
means for better understanding our more sparsely sam¬ 
pled ground-based observations. Examination of the 
EPOXI lightcurve reveals that the shape and ampli¬ 
tude varies from one 17-19 hr cycle (the “single-peaked” 
case) to the next, but that the ensemble of features 
in the lightcurve repeat reasonably well over approxi¬ 
mately 54-55 hr (the “triple-peaked” case). The entire 
EPOXI dataset cannot be satisfactorily phased with a 


Figure 6. Subsect ions of the EPOXI clear filter lightcurve 
(|Williams et al.|2013|) phased to an apparent triple-peaked period. 
The data are identical to those in Figure[5]but are now plotted with 
one color per two days. The colors approximately follow a color 
wheel - green, yellow, orange, red, purple, blue - from October 3 
to November 22. The date range is given in the top left of each 
panel, and the period to which the data are phased is given in the 
top right. Zero phase is set to perihelion (October 28.257). Each 
panel plots exactly two weeks of data with two days of overlap at 
the beginning/end with the panel above/below it. The dates of 
each panel in days relative to perihelion (AT) are —25.0 to —11.0 
(top panel), —13.0 to +1.0 (second panel), —1.0 to +13.0 (third 
panel), and +11.0 to +25.0 (bottom panel). Note that the data be¬ 
come sparse and the amplitude increases during the EPOXI close 
approach (November 4.587; AT = +7.3). 

single period, but subsets of the data can be phased rea¬ 
sonably well and reveal a lengthening apparent period: 
54.6T0.3 hr from October 3-17, 54.8T0.5 hr from Octo¬ 
ber 15-29, 55.0T0.1 hr from October 27-November 10, 
and 55.2T0.2 hr from November 8-22. This is shown in 
Figure [6] where we have normalized the EPOXI data as 
in Figure [5] (including the scaling in order to keep the 
axes the same on both plots), then phased four segments 
to their best by-eye fits with the periods shown on the 
figure and zero phase set to perihelion. We have colored 
the points in two day steps that approximately follow 
the color wheel, with green in early October followed by 
yellow, orange, red, purple, then blue in late November; 
the same colors are used for the same dates in all panels. 

The striking feature of this phased plot is that the 
triple-peaked lightcurve shape does not remain constant 
even during segments that can be phased reasonably well. 
While some features repeat at similar phase and ampli- 
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tude, others evolve steadily from one triple-peaked cy¬ 
cle to the next. Note that this evolution is not due to 
a slight mismatch in the period used for phasing, since 
it occurs in only one area of the lightcurve while other 
areas continue to align well. This can be seen best in 
the bottom panel of Figure [6j showing November 8-22 
(days from perihelion, AT, = +11 to +25). While the 
maxima at ^0.59 and ~0.96 phase remain relatively con¬ 
stant in their shape, height, and location in phase, the 
other maximum changes steadily. The earliest data, in 
light pink/lavender, have the highest peak and reach it 
at a phase of ^0.22. On the subsequent triple-peaked 
cycle (purple), the peak still occurs at a phase of ~0.22, 
but falls dramatically without a significant accompany¬ 
ing change in the height or location of the neighboring 
minima. On the next two cycles (dark blue then royal 
blue/light blue), the maximum at ^0.22 phase contin¬ 
ues to fall and is accompanied by a second maximum 
at ^0.33 phase that increases in height each cycle. The 
final cycle plotted (turquoise) shows only the peak at 
^0.33 phase with no evidence of the former peak at ^0.22 
phase. There are hints that the opposite transition be¬ 
tween the peaks at ^0.33 and ^0.22 phase may be occur- 
ing at the end of the second panel and the beginning of 
the third panel (shades of orange), but gaps in the data 
prevent us from drawing firm conclusions about this. 

The third panel, showing October 27 to November 10 
(AT = — 1 to +13) shows a relatively small amount of 
variation in the lightcurve shape. There is clear evidence 
that the largest measured amplitudes coincide with the 
EPOXI encounter (November 4.587), and progressively 
lessen in the days before and after. The only significant 
variations in the shape are thought to be directly associ¬ 
ated with this. 

The first two panels, showing October 3-17 (AT = 
—25 to —11) and October 15-29 (AT = —13 to +1), 
have sparser coverage than the third and fourth panels, 
but show much more rapid evolution of features in the 
lightcurve. This is most evident at ~0.35 phase in the top 
panel and ^0.62 phase in the second panel, where strong 
peaks are seen just a few triple-peaked cycles after the 
lightcurve was nearly flat at the same phase. The relative 
lack of coherence in the phased lightcurves in these two 
panels is not due to poor choice of the apparent period 
when phasing; these were our optimal apparent periods 
and we could find no apparent period that resulted in 
features aligning as well as they do in the third or fourth 
panels. 

While not shown, we also examined the uncalibrated 
EPOXI data from Septembe r 5-25 provid ed in the sup¬ 
plemental online material of Belton et al. (2013). These 
data were not archived with the PDS because “a small 
light leak at large solar elongations allowed sunlight 
to enter the instrument, causing the comet to appea r 
anomalously bright in all filters” (Williams et al.|[2013). 
In addition, the data were collected much less frequently 
than in October and November. We conducted similar 
analyses to those just discussed for the PDS data, finding 
a shorter apparent period, 54.0+0.5 hr from September 
13-25 (the early September data could not be phased 
conclusively), and more rapid evolution of the triple- 
peaked lightcurve. The much lower sampling rate cou¬ 
pled with the light leak yielded poorer quality results and 
we therefore elected to exclude the September data from 


further discussion. 


5. DISCUSSION 


Having introduced Hartley 2’s coma morphology and 
lightcurve separately, we now consider them in combi¬ 
nation. An obvious question is whether the apparent 
periods determined in each manner are the same, e.g., 
was similar morphology exhibited at similar lightcurve 
phase? We looked at the coma morphology during all 
maxima and minima in our data, since these are the as¬ 
pects of the lightcurve that can be reliably determined 
each cycle despite the lightcurve’s evolution. Within the 
level of uncertainties (<0.5 hr), the same features in the 
coma were seen at the same rotational phase throughout 
the subsets of data that exhibited similar overall coma 
morphology: August through mid-September, early- to 
mid-October, early November, and mid-December. The 
coma’s appearance as a function of position in the 
lightcurve mirrors the the overall morphological evolu¬ 
tion. From August through mid-September the coma 
looked fairly similar at the corresponding lightcurve ex¬ 
trema during all three single-peaked cycles within a 
triple-peaked cycle. In later months, the coma differed 
from cycle to cycle and looked most similar when inte¬ 
ger multiples of three cycles apart. The CN lightcurves 
overlaid with coma morphology for early September and 
early November are shown in Figure [7| 

The correlation between lightcurve phase and morphol¬ 
ogy implies that both techniques should yield compara¬ 
ble apparent periods and that differences are likely due 
to the limitations of the data rather than systematic dif¬ 
ferences due to the method of measurement. This cor¬ 
respondence was not necessarily expected, since projec¬ 
tion effects could cause coma morphology to repeat ear¬ 
lier or later on subsequent cycles. As can be seen in 
Figure [2j previous results from smaller subsets of data 
generally yielded systematically longer apparent periods 
from morphology but, as discussed below, these can now 
be interpreted as due to biases in the limited ground- 
based sampling. Note, however, that this conclusion is 


not relevant for radar observations (e.g., Harmon et al 


2011) as they detected the rotation of the nucleus itself 


rather than variations in the activity or morphology. 

One benefit of plotting the EPOXI data phased to a 
triple-peaked solution is that it revealed that consec¬ 
utive single-peaked cycles were not identical and that 
some portions of the lightcurve evolved over just a few 
days. This was not obvious when plotting single-peaked 
lightcurves from the ground and explains the sometimes 
incompatible portions of the lightcurve we found when 
phasing our more sparsely sampled data. For example, 
the second CN panel in Figure [4] shows phased data from 
early- and mid-September with a strong peak near 0.9 
phase on September 9 and 11 (one triple-peaked cycle 
apart), but no obvious peak at that phase on Septem¬ 
ber 12 (one single-peaked cycle later than September 11) 
or September 2 (10 single-peaked cycles, e.g., one single- 
peaked cycle plus three triple-peaked cycles, earlier than 
September 9). Since the coma morphology was similar 
in all cycles, we were able to determine a robust appar¬ 
ent morphological period, but the apparent period deter¬ 
mined from the lightcurve had much larger uncertainty. 

Phasing the EPOXI lightcurve to a triple-peaked pe¬ 
riod also revealed unequal spacing between lightcurve 
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Figure 7. Triple-peaked CN lightcurve and CN coma morphology for September 1-3 and 9-13 (top) and November 2-8 (bottom). 
The lightcurve details are given in Figure [J] and image details are given in Figure [l] The color schemes are blue/white for Lowell and 
green/white for KPNO, with white representing bright regions and blue/green representing dark regions. The lightcurves were phased to 
the approximate best apparent triple-peaked period at the midpoint of our data. The phase for each image is noted by an arrow color 
coded to match that night’s lightcurve. In both plots the images occurring near lightcurve maxima are shown on the top row and images 
occurring near the lightcurve minima are shown on the bottom row. This reveals that the morphology was much more constant through 
the entire apparent triple-peaked cycle in September than in November. An extremely bright star contaminated all images on September 
9, making the western half of the image appear white and suppressing the coma morphology (although hints of the arc feature can still 
be seen). Also note that September 12 and the second half of September 13 were cloudy; as discussed in the text we have attempted to 
salvage these lightcurve data with comparison stars, but the coma morphology suggests that the lightcurves on these nights are unreliable. 
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peaks. This may contribute to the large scatter in pub¬ 
lished periods since derived single-peaked periods can 
vary significantly depending on which lightcurve peaks 
are used as well as the timing of observations. As an 
example, using the earliest data in the bottom panel of 
Figure[6 (light pink), peaks are seen at ^0.22, ^0.59, and 
^0.96 phase. This yields single-peaked periods of 20.4 
hr, 20.4 hr, and 14.4 hr, all far from the combined mean 
single-peaked period of 18.4 hr. Just such a misleading 
period would be determined if our October 16 and 17 
lightcurves (shown in the third row of the left column 
of Figure [3]) were considered in isolation since the time 
between troughs yields a period of ~21.5 hr (note that 
the coma morphology did not suggest a repetition due 
to the variations in morphology between sin gle-peaked 
cycles). In comparison, Knight & Schleicher (2011) de¬ 
termined an apparent period of 18.15iU.15 hr around 
this time based o n coma morphology over several nights, 
while |Belton et al. (2013) determined an apparent period 
of 18.02±0.06 hr from the EPOXI lightcurve. 

Shifts in the location in phase space of a lightcurve 
feature over a short interval, as seen in the EPOXI 
lightcurve when phased to the triple-peaked solution, can 
also distort the inferred single-peaked period. Returning 
to the bottom panel of Figure [6j periods determined us¬ 
ing the first peak in purple (at ~0.22 phase) and royal 
blue (at ^0.33 phase) just two triple-peaked cycles apart 
would imply a single-peaked period of 19.4 hr for the 
six intervening cycles, significantly longer than the com¬ 
bined mean single-peaked period of 18.4 hr. The EPOXI 
lightcurve reveals a number of such transitions, and oth¬ 
ers undoubtedly occurred during the gaps in EPOXI cov¬ 
erage and before/after EPOXI observations. 

In the two preceding cases, averaging over a large num¬ 
ber of intervening cycles will reduce the scatter in the in¬ 
ferred single-peaked period. This explains the relatively 
smooth progression of single-peaked periods across the 
apparition that we determined from coma morphology 
on consecutive observing runs since these were separated 
by tens of intervening single-peaked cycles. However, 
single-peaked periods measured from just a few days of 
data can be heavily biased. This may explain much of 
the scatter in published solutions, since most were de¬ 
rived from data taken over just a few days. 

A related uncertainty in the inferred single-peaked pe¬ 
riod comes from imperfect sampling from the ground due 
to the diurnal cycle. We typically observed the comet for 
5-9 hr per night over 3-5 consecutive nights of a single 
observing run. Due simply to the cadence of our ob¬ 
servations and the separation between observing runs, 
we often got little to no data on one of the three cycles 
that make up the triple-peaked period. We examined 
this by calling the three single-peaked cycles that make 
up one three-peaked cycle A, B, and C, arbitrarily start¬ 
ing A at 0:00 UT on the first night of an observing run, 
and defining it as extending for one single-peaked pe¬ 
riod (using the period we previously determined for that 
epoch), at which point cycle B started, followed by C, 
then by A, etcH While we had decent sampling within 
all three cycles tor our first two pairs of observing runs 

7 These cycle start/stop times are set by the da tes of our ob¬ 
servat ions and were not the same as those used by |Drahus et al.| 
([201 lj to partition their data. 


(mid-August to early September and early September 
to mid-September), after mid-September cycle C was al¬ 
most never imaged. As just discussed, this effect is mit¬ 
igated by the large number of intervening cycles, but 
still might yield inferred periods biased by ^0.1-0.3 hr 
in widely spaced data. It is difficult to quantify how 
this affects our paired December images since EPOXI 
was no longer observing so we do not have full coverage 
of the triple-peaked lightcurve, but given that our pairs 
were separated by five or fewer single-peaked cycles, the 
uncertainty could be 1 hr or more. Note that the un¬ 
certainties listed in Table [2] are the standard deviation 
of our measured periods and do not account for such bi¬ 
ases, nor do the range of viable apparent periods given 
in Table |31 

We tested if the unequal spacing between lightcurve 
peaks observed in EPOXI data from October and Novem¬ 
ber existed earlier in the apparition as follows. We took 
the interval between a matched pair of images and sub¬ 
tracted the apparent period determined for that epoch 
times (TV —1) where N is the number of intervening pe¬ 
riods we determined for that image pair in Section [3J If 
the length of the single-peaked cycles was uniform, this 
“remainder” should be approximately equal to the ap¬ 
parent period; if unequal the remainder will be longer 
or shorter depending on the starting and ending single- 
peaked cycles. We found remainders that differed from 
the apparent period by increasing amounts during the 
apparition, from about 1 hr in the mid-August to early 
September image pairs to as much as 7 hr in mid-October 
to early November image pairs. The large range of re¬ 
mainders should not be interpreted as being due to large 
uncertainties in the matching of coma morphology, as we 
see clustering in the remainders having the same start¬ 
ing and ending single-peaked cycles. That is, remainders 
determined from an initial image in cycle A (as defined 
in the previous paragraph) and a final image in cycle B 
were similar, remainders determined from an initial im¬ 
age in cycle C and a final image in cycle A were similar 
to each other but were different from the A-B pair re¬ 
mainders, etc. While all of the biases discussed above 
may skew the specific remainder values, the increasing 
range of remainders strongly suggests that the length of 
the single-peaked cycles were more equal in August and 
early September than in October and November. This 
finding extends earlier than the EPOXI monitoring and 
therefore is a valuable constraint for future efforts to un¬ 
derstand the component periods of NPA rotation (which 
we presume are responsible for the observed behavior). 

Despite the large uncertainties in our periods after 
early November, our combined data suggest the apparent 
period chan ged little afte r the E POXI flyby. This is con¬ 
sistent with Belton et al.| ([2013])’s analysis of the EPOXI 
spacecraft data which founcTa period of 18.63T0.02 hr 
on November 11 and 18.61T0.09 hr on November 20. In 
light of our better understanding of how sampling bias 
can skew period determinations, it is also probably con¬ 
sistent with our earlier estimate of an apparent period 
“near 19 hr” on December 9-10 (Knight & Schleicher 
2011) and|Jehin et al. (|2010 )’s report that the apparent 
period increased from 18.2 to 19 hr from October 29 to 
December 7 (the 0.3 hr error bar on their period from the 
first two weeks of data, when the comet was brighter, im¬ 
plies their uncertainties by early December were at least 
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this large). 

The gas production rates of Hartley 2 pe aked 0-20 days 
post-perihelion (e.g., Knig ht fe Schleicher| 2013). Assum¬ 
ing a simple model where the outgassing torq ue is pro¬ 


por tional to the gas production rate, Knight fc Schle- 
icher (2011) argued for non-linear changes in the appar- 


ent period, and a decrease in the rate of change of the 
period was predicted starting from early November (see 
their Figure 6). However, the extent of the flattening 
of the rate of change of the period appears to be much 
more than suggested there. The rate of change of the 
rotational state, including those of the component peri¬ 
ods, depends on the net component torques acting on the 
nucleus due to outgassing. Torques can be due to out- 
gassing from specific source regions on the surface or due 
to much more widely distributed activity on the surface 
of the nucleus. Regardless of the primary cause, the out¬ 
gassing pattern and the resultant net torques can easily 
change due to different insolation patterns as a comet 
moves along in its orbit. In other words, the net com¬ 
ponent torques acting on the nucleus are not necessarily 
proportional to the gas prod uctio n rate or to th e level 
of surface activity (cf. Samarasinha & Mueller 2013). 
Due to this complexity, one should be cautious in di¬ 
rectly correlating this flattening observed for the period 
in November and December to the activity of the nu¬ 
cleus. The apparent disconnection between the torquing 
and the production rates is another constraint for future 
modeling efforts. 

We found a steady (within the uncertainties) increase 
in the apparent triple-peaked period when phasing the 


EPOXI data , in c ontrast to ear l ier an alyses by | A’Hearn 
et al.| (|2011|) and |Belton et al.| (|2013|) who found more 
if 


variability in the period during the same interval. This 
is likely due to how each analysis technique handles 
the most rapidly evolving segments of the lightcurve. 
When visually assessing goodness of fit, we were cog¬ 
nizant of the evolution of the lightcurve from cycle to 
cycle. Thus, we disregarded regions in obvious transi¬ 
tion, instead focusing on those regions that remained 
relatively unchanged during a given interval. Since the 
lightcurve was continuously evolving, focusing on align¬ 
ment of features in the more rapidly evolving regions 
would have yielded much larger uncertainty and, likely, 
a more chaotic trend in the apparent period during the 
EPOXI observing window. In contrast, the power spec¬ 
tra technique em ployed by |A’Hearn et al .| (|2Q11|) and 
Belton et al. (2013) treats all data equally, pelton et al. 
(|2U13|) found very small uncertainties (<0.07 hr) when 
the lightcurve was most stable, in late October and early 
November, but found apparent periods with more varia¬ 
tion and considerably larger uncertainties (0.86-5.76 hr) 
in September and early October, when evolution of the 
lightcurve was significant, the data were sampled less fre¬ 
quently, and (in September) there was a light leak that 
varied over time. The large uncertainty in September 
makes it difficult to determine a clear trend in this com¬ 


pone nt peri od in either A’Hearn et al. (2011 )’s or Belton 
et al. (2013)’s analyses. We surmise that Hartley 2’s com- 


ponent periods changed steadily during the apparition 
and therefore conclude that the apparent triple-peaked 
periods we determined from visually phasing the data are 
more plausible. 

The fact that the triple-peaked and single-peaked ap¬ 


parent periods are almost exactly in a 3:1 ratio through¬ 
out the EPOXI observing window should not be surpris¬ 
ing - the triple-peaked apparent period has three peaks 
and by definition should yield a single-peaked apparent 
period that is 1/3 as long. What is surprising is that 
two of the underlying component periods of NPA rota¬ 
tion appear to have been close enough to a 3:1 ratio (or 
some alias such as 3:2) that the triple-peaked lightcurve 
exhibited minimal evolution during the third panel of 
Figure |6l (October 27 to November 10). The lightcurve 
clearly showed more rapid evolution in the surrounding 
weeks (shown in the first, third, and fourth panels of 
Figure pi). This seems to imply that at least one of the 
component periods of NPA rotation was changing such 
that their relative ratio was approximately an integer 
multiple sometime between October 27 and November 
10, and was further away from an integer multiple both 
earlier and later. A simil ar con clusion is reached using 


the ratio of Belton et al. (2013)’s and zq, which were 
in a 3:1 ratio approximately October 27. The transition 
in lightcurve shape is a hallmark of NPA rotation, and 
similar behavior was se en in the lightcurve of lP/Halley 
(Schleicher et al.| 1990 Schleicher et al., in prep.). This 
serves to - remind us that the apparent periods we are 
measuring are not necessarily the component periods of 
the NPA rotation, but the interaction of the underly¬ 
ing periods manifesting themselves as repetition in coma 
morphology and brightness. 

An additional reminder of the evolving ratio of the 
component periods is the variation in morphology at 
the the same lightcurve extrema seen in Figure [7| In 
September, the morphology near all three peaks in the 
lightcurve is similar, as is the morphology near all three 
troughs. In contrast, the morphology varies from one 
peak to the next in November and the bulk brightness 
flips hemispheres during the middle cycle. While beyond 
the scope of the current paper, this behavior will be an 
important constraint for future modeling efforts. 

6. SUMMARY 

In this work we have combined observations of Hart¬ 
ley 2 during the 2010 apparition obtained at Kitt Peak 
National Observatory and Lowell Observatory. The 
two datasets were not contemporaneous except at the 
time of the EPOXI flyby, and the combined dataset 
roughly halves the time between consecutive observing 
runs as compared to data collected at a single observa¬ 
tory. This allowed determinations of intermediate pe¬ 
riods that were not possible between consecutive runs 
with the same telescope. We measured the apparent pe¬ 
riod of CN coma morphology repeatability between runs, 
yielding apparent periods additional to those previousl y 
reported from our individual datase ts (Knight & Schle- 
icher||2011| |Samarasinha et al.||2011 ). VVe also measured 
inner coma lightcurves in both continuum (R-band) and 
gas (CN) within and between observing runs, the first 
such lightcurves of Hartley 2 obtained from the ground 
to be published. Phasing of these lightcurves yielded 
additional apparent periods and revealed that the coma 
morphology and lightcurves were in sync. 

Our new periods are consistent with the roughly linear 
increase in apparent period reported by previous authors 
up to mid-November. The post -EPOXI encounter peri¬ 
ods suggest that the period flattened out or possibly even 
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decreased slightly in late November and December. The 
next favorable observing window for Hartley 2 will occur 
in mid-2016, prior to the onset of strong activity. Ob¬ 
servations at this time should reveal the current period 
and test how it has changed since late 2010. However, 
as this will still be at relatively large heliocentric dis¬ 
tances, the brightness may be dominated by the nucleus, 
in which case large aperture telescopes will be required 
and the lightcurve may display different characteristics 
than were seen in the activity dominated lightcurve of 
2010. Unfortunately, Hartley 2 next reaches perihelion 
on the far side of the Sun in April 2017, so Earth-based 
observations designed to look for further changes in the 
period will likely not be possible until early 2018 and will 
again require large aperture telescopes. 

We also examined the EPOXI lightcurve in order to 
better understand our own data. Phasing of the EPOXI 
data to a triple-peaked period revealed that consecu¬ 
tive single-peaked cycles had different lightcurve shapes 
and separations between peaks that approximately re¬ 
peated every third cycle. Segments of the triple-peaked 
lightcurve evolved from one triple-peaked cycle to the 
next while other segments remained approximately con¬ 
stant for many cycles. These behaviors can explain much 
of the scatter between apparent periods measured at sim¬ 
ilar epochs in our own and other datasets. We showed 
that limitations of undersampled datasets can lead to dif¬ 
ferences from the true period of 1 hr or more, making it 
more important to have longer baselines. 

Despite Hartley 2 having been extensively studied from 
the ground and space in support of EPOXI , its rota¬ 
tion is not yet fully understood. The evolution of the 
lightcurve shape and coma morphology seem to indicate 
that the relative ratios of the component periods of NPA 
rotation changed during the apparition. However, spe¬ 
cific change(s) corresponding to component periods are 
currently not known. At present there is also no consen¬ 
sus rotational angular momentum vector (RAMV), with 
multiple author s suggesti ng very different RAMYs (e.g., 
A’Hearn et al7||201H [Harmon et al7||2Qll| IWaniak et al 


2012| Knight &; Schleicher|2Q13||Belton et al.|2Q13| ). Fur¬ 

thermore, the type of JNPA rotation (LAM versus SAM) 


(LAYT 


SAM) 


is not yet conclusively known. Hartley 2’s NPA rotation 
is challenging to model due to the large number of free 
parameters, but numerical simulations may yet arrive at 
a solution that matches the coma morphology and repli¬ 
cates the complex lightcurve. Such a result would then 
be very useful for planning a future spacecraft mission to 
this compelling comet. 
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